val_flex <- mean(clean_LISS_data$high_flex, na.rm = TRUE)
val_flex_aftcov <- mean(clean_LISS_data$high_flex_aftcov, na.rm = TRUE)

val_tel <- mean(clean_LISS_data$high_telecommute, na.rm = TRUE)
val_tel_aftcov <- mean(clean_LISS_data$high_telecommute_aftcov, na.rm = TRUE)

covid_table <- data.frame(col1 = c(val_flex, val_flex_aftcov),
                          col2 = c(val_tel, val_tel_aftcov))

rownames(covid_table) <- c("Before Pandemic", "After Pandemic")
colnames(covid_table) <- c("High Flexibility", "Telecommuting")

table_covid <- xtable::xtable(covid_table)
writeLines(xtable::print.xtable(table_covid, include.rownames = TRUE, include.colnames = TRUE, floating = FALSE, type = "latex", sanitize.colnames.function = identity), file.path(graph_dir, "table_covid.txt"))